Effect of inelastic collisions on multiphonon Raman scattering in graphene 
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We calculate the probabilities of two- and four-phonon Raman scattering in graphene and show 
how the relative intensities of the overtone peaks encode information about relative rates of different 
inelastic processes electrons are subject to. If the most important processes are electron-phonon and 
electron-electron scattering, the rate of the latter can be deduced from the Raman spectra. 



Introduction. — Collisions of Dirac electrons are quali- 
tatively different from those of electrons in conventional 
metals: energy and momentum conservation leave no 
phase space for relaxation of a quasiparticle excited above 
the Dirac vacuum. Electron collisions in graphene and re- 
lated compounds continue to be a subject of theoretical 
studies 0, [^jlli]. Thus, any experimental information on 
collisions of Dirac electrons would be extremely valuable. 
However, to experimentally separate contributions from 
different mechanisms to electron lifetime is a hard task. 
The present work suggests a way to separate electron- 
phonon and electron-electron contributions to the elec- 
tron lifetime by analyzing Raman spectra. 

Raman spectrum of graphene consists of distinct peaks 
corresponding to different optical phonon branches as 
well as their overtones. Thus, Raman scattering measure- 
ments represent a powerful experimental tool for study- 
ing phonon modes, as well as their interaction with elec- 
trons (since electronic excitations are involved in Ra- 
man scattering as intermediate states). Indeed, electron- 
phonon interaction and Raman scattering in graphene 
has attracted a great deal of interest, both experimen- 
tal 0, H, S 0, il and theoretical [i, 0^. Here we show 
how even more information can be extracted from Raman 
spectra. 

Qualitative picture. — Photon wave vector is negligible, 
so momentum conservation requires that Raman scatter- 
ing on one intervalley phonon must be impurity-assisted 
[process (b) in Fig. [U giving rise to the so-called D Ra- 
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FIG. 1: Schematic representation of the role of electron dis- 
persion (Dirac cones, shown by solid lines) in Raman scat- 
tering: (a) one-phonon G peak, (b) impurity-assisted one- 
phonon D peak, (c) two-phonon D* peak. Vertical solid ar- 
rows represent interband electronic transitions accompanied 
by photon absorption or emission (photon wave vector is ne- 
glected), dashed arrows represent phonon emission, the hori- 
zontal dotted arrow represents impurity scattering. 



man peak] . D peak is absent in the experimental Raman 
spectrum of graphene [3| , so impurity scattering is negli- 
gible in these samples, and is disregarded hereafter. 

Looking at the intermediate electronic states (Fig.[T]), 
we notice that for one-phonon scattering [processes 
(a), (b)] at least one intermediate state must be vir- 
tual, since energy and momentum conservation cannot 
be satisfied simultaneously. For the two-phonon scatter- 
ing [process (c)] all electronic states can be real. We em- 
phasize the qualitative difference between the fully reso- 
nant process (c) and the double-resonant [lH process (b) , 
where one intermediate state is still virtual. 

Obviously, our argument can be extended to all multi- 
phonon processes with odd and even number of phonons 
involved: in order to annihilate radiatively, the electron 
and the hole must have opposite momenta; if the total 
number of emitted phonons is odd, the electron and the 
hole must emit a different number of phonons, which is 
incompatible with energy conservation in all processes. 

Real even-phonon processes can be viewed in the fol- 
lowing way. The incident photon creates an electron 
and a hole - real quasiparticles which can participate 
in various scattering processes. If the electron emits a 
phonon with a momentum q, the hole emits a phonon 
with the momentum — q, and after that the electron and 
the hole recombine radiatively, the resulting photon will 
contribute to the two-phonon Raman peak. If they do 
not recombine at this stage, but each of them emits one 
more phonon, and they recombine afterwards, the result- 
ing photon will contribute to the four-phonon peak, etc. 

Besides phonon emission and radiative recombination, 
electron and hole are subject to other inelastic scatter- 
ing processes, which can also be viewed as emission of 
some excitations of the system. The key point is that 
for real quasiparticles, the probability to undergo a scat- 
tering process a is determined by the ratio of corre- 
sponding scattering rate 2^^ to the total scattering rate 
27 = X)q27c(, not by history. This probability deter- 
mines the relative frequency-integrated intensity of the 
corresponding feature in the Raman spectrum. Thus, 
the ratio of integrated intensity /2n+2 of the Raman peak 
corresponding to 2n -I- 2 phonons to that for 2n phonons 
(/2„) must be proportional to {'^ph/lY, where 2^ph is 
the rate of emission of each of the two phonons, and the 
square comes from the phonon emission by the electron 
and the hole. 
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In the Raman spectrum of graphene two two-phonon 
peaks are seen: the so-called D* peak near 2u!Ai = 
2650 cm^^, and the G* peak near 2uje2 = 3250 cm^^, 
corresponding to scalar Ai phonons from the vicinity of 
the K point of the first Brillouin zone, and to pseudovec- 
tor E2 phonons from the vicinity of the F point, respec- 
tively. The D* peak is more intense, thus the most inter- 
esting for practical purposes is to compare the intensities 
of D* and its four-phonon overtone at 5300 cm~^, which 
we win call 2D*: 



ho'/lD' =0.14(7A,/7) 



(1) 



The coefficient 0.14 was obtained by direct calculation 
assuming Wi„ ^ ujai ^ 7i where uji^ is the incident 
photon frequency. 2jai is the rate of emission of the 
Ai phonon. It depends on the electron energy, which can 
be taken to be ujin/2. Eq. ([T]) represents the main result 
of this paper. Let us now discuss its practical implication. 

In graphene, the most obvious competitor of the 
phonon emission is the electron-electron scattering: the 
optically excited electron can kick out another one from 
the Fermi sea, i. e., to emit another electron-hole pair 
(for an electron above the Dirac vacuum the phase space 
for an intravalley collision is zero, so the collision has 
to be intervalley or impurity-assisted [3). Thus, Raman 
spectrum should contain contribution from electron-hole 
pairs; however, their spectrum extends all the way to the 
energy of the photo-excited electron (optical energy) in 
a completely featureless way. Thus, it cannot be distin- 
guished from the parasitic background which is always 
subtracted in the analysis of Raman spectra, and cannot 
be seen in the Raman spectrum directly. However, as- 
suming 7 = + IE2 + 7ee, where 2'~^e2 is the rate of 
the E2 phonon emission, and 27ee is the electron-electron 
collision rate, one can extract the value of 7ee from the 
experimental data using Eq. ([1]) , relative to phonon emis- 
sion rates. More precisely, in this way one obtains the 
rate of all inelastic scattering processes where the elec- 
tron loses energy far exceeding the phonon energy. 

Note that arguments leading to l2n+2/l2n oc (riph/lY 
are not specific for graphene; in fact, this is nothing but 
Breit-Wigner formula, applied once for the electron and 
once for the hole. Multi-phonon Raman scattering has 
been studied in wide-gap semiconductors both experi- 
mentally 12, [1^ (up to ten phonons were seen in the 
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Raman spectra of CdS), and theoretically 
a wide-gap semiconductor an optically excited electron 
does not have a sufficient energy to excite another elec- 
tron across the gap, so the electron-electron channel is 
absent. In addition, interaction with only one phonon 
mode is dominant, so the ratios of subsequent peaks are 
represented by a sequence of fixed numbers. The sim- 
ple band structure allowed a calculation of the whole se- 
quence. For graphene, we restrict ourselves to the calcu- 
lation leading to Eq. 
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TABLE I: Irreducible representations of the group Csv and 
their characters. 
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valley-off-diagonal matrices 
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TABLE II: Classification of 4x4 hermitian matrices according 
to irreducible representations of the Ce^ group. 



Model. — We measure the single-electron energies from 
the Fermi level of the undoped (half- filled) graphene. The 
Fermi surface of undoped graphene consists of two points, 
called K and K' . Graphene unit cell contains two atoms, 
each of them has one 7r-orbital, so there are two electronic 
states for each point of the first Brillouin zone (we dis- 
regard the electron spin). Thus, there are exactly four 
electronic states with zero energy. An arbitrary linear 
combination of them is represented by a 4-component 
column vector if). States with low energy are obtained by 
including a smooth position dependence '4'{v), r = {x,y). 
The low-energy hamiltonian has the Dirac form: 



H^i^ j Sri;\v) (-i-yS- V)V^(r). 



(2) 



We prefer not to give the explicit form of the isospin ma- 
trices S = (Ea;, Ey), which depends on the choice of the 
basis (specific arrangement of the components in the col- 
umn ij}). We only note that all 16 generators of the SU (4) 
group, forming the basis in the space of 4 x 4 hermitian 
matrices, can be classified according to the irreducible 
representations of Cg^ , the point group of the graphene 
crystal (Tables |T] and |Tl|. They can be represented as 
products of two mutually commuting algebras of Pauli 
matrices Ea;,E„,E2 and A^ 



,Ay,K^ 



3 12 1, which fixes 



their algebraic relations. By definition, Ex,Ey are the 
matrices, diagonal in the K, K' subspace, and transform- 
ing according to the Ei representation of Cgu . 

We restrict our attention to scalar phonons with wave 
vectors close to K and K' points - those responsible for 
the D* Raman peak. The two real linear combinations 
of the modes at K and K' points transform according 
to Ai and Bi representations of 0%^ and are shown in 
Fig. O We take the magnitude of the carbon atom dis- 
placement as the normal coordinate for each mode, de- 
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FIG. 2: Phonon modes responsible for the D* Raman peak. 



noted by Ua and Uh, respectively. Upon quantization of 
the phonon field, Ua, Ub and the phonon hamiltonian i/ph 
are expressed in terms of the creation and annihilation 
operators 6^^, 6q^, /i = a, 5, as 



h.c. 



y/2NMuJAi 

(3) 

The crystal is assumed to have the area L^Ly, and to 
contain N carbon atoms of mass M. The q summation 
is performed as ~^ ^xLy J d^q/(27r)^. "h.c." stands 
for hermitian conjugate. By symmetry, in the electron- 
phonon interaction hamiltonian [l8| the normal displace- 
ments Uf^ couple to the corresponding valley-off-diagonal 
4x4 matrices from Tabic [TTJ 



H^nt = 3Fa, J d'rV^W [ti„(r)A,S, + Ub(r)Ay^,] ^(r). 

(4) 

Interaction with light is obtained from the Dirac hamil- 
tonian (12) by replacement V ^ V — i{e/c)A, where the 
vector potential A is expressed in terms of creation and 
annihilation operators aQ^,aQ.£ of three-dimensional 
photons in the quantization volume V = L^LyL^, la- 
beled by the wave vector Q and two transverse polariza- 
tions £ ~ 1,2 with unit vectors eq^^: 



Q, 



2ttc 



iQr 



VQ 



h.c.) 



(5) 



Calculation. — A^(qi, /ii; . . . ; q„, /i„), the matrix ele- 
ment of the transition from the state with one incom- 
ing photon (frequency LOi„, polarization ei„) into the 
state with one outgoing photon (frequency uJout, polariza- 
tion Bout) and n phonons (modes /ii, . . . , /i„, wave vectors 
qi, . . . , q„), is calculated as shown in Fig. [3l We empha- 
size the necessity to include the inelastic broadening 7 
in the electronic Green's functions, as the dominant con- 
tribution to the integral comes from regions where the 
denominators are small (~ 7, corresponding to real elec- 
trons and holes, as discussed above). Given the matrix 
element, we sum over final photon and phonon states, 
and express the absolute dimensionless probability of n- 
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FIG. 3: (a), (b) Fully resonant diagrams giving the matrix 
element M for the two- and four-phonon Raman scattering, 
respectively. Other possible diagrams (not shown) are not 
fully resonant, and are smaller by a factor 'y/cuAi ■ (c) Ana- 
lytical expressions, corresponding to lines and vertices of the 
diagrams. In the electron loop integration over the internal 
momentum and energy, as well as tracing over 4x4 matrix 
structure should be performed. 



phonon Raman scattering as 

X E E 

qi + ...+q,.=0{;i,} 



^M(P{q„/.J) 



, (6) 



where V denotes permutations of phonon arguments. 
For the D* peak the diagram in Fig. [3^ gives 



X(q) 



X Ginjz [eq X 6^ 



v{q - qts) - 2ij]^/'^ 

(7) 

t^o«t)/(2i'), 



where Gq = q/|q|. The value — (wj 
around which expression ([7]) is strongly peaked, corre- 
sponds to backscattering of the electron and the hole by 
the phonons. This sharply peaked dependence cannot 
be obtained from pure symmetry considerations, which 
just prescribe vanishing of the matrix element when the 
scattering angle (p = [l^. Its origin is the quasiclassi- 
cal nature of the electron and hole motion [201 ; the dis- 
persion of tp being determined by quantum diffraction: 
— 7r| ~ y^j/ujin <C 1. One consequence of this peaked 
dependence is that the width of the D* peak is deter- 
mined by the electron and phonon lifetimes only, not by 
the phonon dispersion. Besides, it should lead to a signifi- 
cant polarization anisotropy of the D* peak [131 . Here we 
simply sum over the polarizations; substituting Eq. ([T]) 
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FIG. 4: The self-energy diagram describing the phonon emis- 
sion. The dashed hne represents the phonon Green's function 



prescription. Substituted in Eq. ([T]), it gives 'fee ~ 0.447, 
which agrees with the calculated 27ee « 10 meV Q and 
the total 2j measured by time-resolved photoemission 
spectroscopy (20 meV in Ref. 25 meV in Ref. [11], 
all values taken for e = Win/2 = 1 eV). Even though a 
recent ARPES measurement gives a significantly larger 
value for 27 ^ 100 meV [2^, the issue of validity of the 
tight-binding model for electron-phonon coupling seems 
to deserve further investigation. 



into Eq. 



Id- 



and approximating LUout 
(e^/c)^ v'^ Lo. 



487r 



j2 

in 



q2 ry2 



9Fl 



we obtain 
27a2 



(8) 



We can allow for trigonal warping and electron-hole 
asymmetry terms in the dispersion of electrons and holes: 
zLvp ±{vp + Q!3p2 cos 3(/3p) + aop^, where ipp is the 
polar angle of p, tight-binding model gives 0:3 = wa/4 
{v w 10^ cm/s, a « 1.42 A), and ao(l eY)^/v^ - 
0.1 eV [2l[. The relative corrections to expression ([8]) 
from these terms are (9/8)(a3WAi/w^)^ ~ 10"'' and 
-[aQ{u,nLOAjv^)/{2j)]y2 ^ 10-4(c^,„/27)2. 

Evaluation of the diagram in Fig. [3b gives the intensity 
of the four-phonon Raman peak: 



/2£,. = 0.088 



647r^ 7^ 



MloaiV' 



(9) 



Finally, we express the Raman scattering probabilities 
in terms of the phonon emission rate. The latter is given 
by the imaginary part of the self-energy (Fig. 0]): 



1A^ (e) 



27a' 



Eqs. 



(10) 



and Eq. (fTO]) 



where 9{e) is a step function, 
with e = LUin/2 give Eq. ([1]). 

Instead of conclusion, wc quote an experimental value 
for the ratio Id- /ho* ~ 40 [2^1, so Eq. ^ gives 
lAi/l ~ 0.42. We assume 7 = 7^1 + 7B2 + 7ee and 
note that the emission rate 27^2 of E2 phonons is de- 
scribed by Eq. pO|) . with the replacements Fa^ Fe^ 
(the corresponding coupling constant), loa^ — > w^Jj. In 
the tight-binding model Fa^ = Fe2 (for the normaliza- 
tion of the phonon displacements chosen here), which 
agrees with a DFT calculation 23] up to 1%. The as- 
sumption Fai — Fe2 gives 7ee ~ 0.22 7. 

On the other hand, Fa-^ and Fe2 are not related by 
any symmetry. For intensities of the two-phonon Ra- 
man peaks D* and G* our calculation gives Id* /Ic — 
2(Fa, /Fp^^YiijjEp/i jJA, Y ■ the experimental value being 
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Id'/Ig' « 20 [4,122]. This suggests FaJFe^ « 1-6, 
in significant disagreement with the tight-binding model 
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